Talk given at the Workshop "RADIATIVE CORRECTIONS relevant for the 
HERMES EXPERIMENT", Zeuthen, Berlin, Germany, 21-23 June, 1993 



Radiative Corrections for HI F 2 (x,Q 2 ) Measurement 



A. A. Akhundov 
DESY - Institut fur Hochenergiephysik Zeuthen, 
Platanenallee 6, D- 15738 Zeuthen, Germany 
and 

Institute of Physics, Azerbaijan Academy of Sciences 
pr. Azizbekova 33, 370143 Baku, Azerbaijan 

ABSTRACT 

The numerical study of QED radiative corrections in the low x region for the 
HI experiment at HERA was performed. The analytical program TERAD91 
was used to get the size of the radiative corrections for the measurement of 
the proton structure function F 2 {x,Q 2 ) in the range x = 10 -2 — 10 -4 and 
5 GeV 2 < Q 2 < 80 GeV 2 . It was found that radiative corrections can be not 
small in some points of the phase space of the mixed variables. 

In the long period from the SLAC experiments 1 till the operation of HERA, the 
experiments with neutral current deep inelastic ep-scattering 

e(h)+p( Pl ) ->e(k 2 )+X(p 2 ) (1) 

had a relatively simple structure, and the measurements completely rested on the 
registration of energy and angle of the scattered electron. The new detector gener- 
ation, the HERA detectors HI 2 and ZEUS, 3 allows to measure both the scattered 
electron and hadronic final state. 

For the physical analysis of the ep-collisions at HERA one can use not only the 
familiar electron variables 

Ql = -(h - h) 2 , y e =Pi(k 1 -k 2 )/p 1 k 1 , x e = Q 2 e /(sy e ), (2) 

where 

s = {h + pi) 2 = 4E e E p , (3) 
but also the kinematical variables from the hadron measurement 



Ql = ~(P2 ~ Pif, Vh = Pi(P2 - Pi)/pih, x h = Ql/ (sy h ), (4) 



or some mixture of both. Here E e and E p are the energies of incident electron and 
proton. 

In order to extract from experimental data the inclusive cross section of deep 
inelastic ep-scattering in one-boson-exchange approximation one must to take into ac- 
count contributions from higher order QED 4 and electroweak processes. For example 
the bremsstrahlung process: 

e(h) + p(pi) - e(k 2 ) + X(p 2 ) + 7 (fc), (5) 

and also the elastic radiative tail : 

e(fci) + p(pi) -> e(fc 2 ) + p(p2) + 7(*0- (6) 

contribute to the observed cross section of the reaction ( 1). 

For different choices of the variables there are differences in the predictions for 
the contribution of the higher order processes. 

The theoretical investigations of these radiative effects, or radiative corrections 
(RC), for experiments at HERA were summarized in Proceedings of the Workshop 
"Physics at HERA" . 5 Different programs to calculate the radiative corrections have 
been cross-checked . 6 

The analysis of neutral current deep inelastic scattering data in the HI experi- 
ment at HERA in 1992 7 was based on two independent approaches: 

I) " detector oriented" approach ( " Brussels- Paris-Saclay analysis" ) 

II) "physics oriented" approach ( "Zeuthen analysis" ) 

Within the analysis I the kinematics of the inclusive deep inelastic scattering 
process (1) were determined from the energy E' e of the scattered electron, and its 
polar angle 9 e , measured relative to the proton beam direction 

Q 2 e =4E e E' e cos 2 (9 e /2), (7) 



y e = l-(E' e /E e )sm 2 (9 e /2). (8) 

The cross sections, acceptances of detector, radiative corrections and resolution 
effects were obtained in (^J~E~l,6 e ) bins. The cross sections, that have been measured 
using this binning, were transformed into cross sections in (x e , Ql)-bms. 

In the analysis II the kinematics of the process (1) were obtained from a 
measurement of the electron variables E' e , 9 e and scaling variable yn from the hadrons. 
The hadronic variable (4) can be determined using the Jacquet-Blondel method 8 



(9) 



where is the energy of a hadron and p z ^ its momentum component along the 
incident proton direction. 

The momentum transfer Q 2 was determined from the electron variables 9 e and 
E' e as in (7) and Bjorken x by using a mixture of the electron and hadron measure- 
ments 9 

Xm = Q 2 e /(sy h ). (10) 

In the analysis II cross sections, radiative corrections and efficiencies were cal- 
culated in (x m , Qe)-bins. 

A further difference between the two methods of the analysis of the data was 
the use of two different types of programs — the Monte Carlo event generators HER- 
ACLES 10 in the analysis I and the analytical program TERAD91 11 in the analysis 
II — for the calculations of the radiative correction factor S(x,Q 2 ). The radiative 
correction factor 5(x, Q 2 ) is defined by the equation 

2 _ d 2 a theor /dxdQ 2 
[X,Q } ~ d 2 a B °™/dxdQ 2 ~ 

where d 2 a theor /dxdQ 2 is the theoretical approximation to the measured cross section 
d 2 a meSLS / dxdQ 2 , which contains the contributions from higher order electroweak pro- 
cesses, and d 2 a Born /dxdQ 2 is the Born cross section of the process (1). 

In this note we present the numerical calculations of the RC factor 5(x, Q 2 ) 
for the "Zeuthen analysis" of the proton structure function F 2 (x, Q 2 ) in the range 
x = 10- 2 - 1(T 4 and 5 GeV 2 < Q 2 < 80 GeV 2 in the HI experiment at HERA 7 

The Monte Carlo HERACLES is suited for the calculation of RC to the differ- 
ential cross sections in terms of the electron variables because the program calculates 
cross sections and generates events in terms of the variables x e and Q 2 . Therefore 
most of the generated radiative events in HERACLES have low Q\ and the calcula- 
tion of the radiative corrections in the terms of other variables ( hadron, mixed ) can 
be performed with less efficiency only. For this reason in the analysis II was used the 
analytical program TERAD91. 

TERAD91 - a program package to calculate radiative corrections for deep in- 
elastic ep-scattering - was created in 1991 during Workshop "Physics at HERA" and 
accumulates about 15 years experience in the field of calculation of the radiative cor- 
rections for deep inelastic ZiV-scattering . 12 The program TERAD91 can be applied to 
calculate the RC factor of the order 0(a) to the measured differential cross section, 
d 2 a meas /dxdQ 2 , after the acceptance and smearing corrections. 

The program package TERAD91 contains one part, TERAD, which is based 
on the model-independent treatment 13 of the leptonic QED corrections for neutral 
current ep-scattering, including bremsstrahlung process ( 5) and QED vertex cor- 
rections. TERAD does not necessarily rely on the quark parton model and uses the 
phenomenological structure functions F 1 ,F 2 and F 3 , which are defined in the Born 
approximation. 



(11) 



The program TERAD91 can be used for the model-independent calculations of 
QED corrections for neutral current deep inelastic ep-scattering in the terms of several 
variables - electron, hadron, mixed variables 14 - and in the new version TERAD93 15 
also in terms of the Jacquet-Blondel variables . 16 The leptonic QED corrections are 
dominated and must be calculated carefully. 

In the low Q 2 region the Born cross section of the deep inelastic scattering (1) 
is determined by two structure function F 2 and 2xFi = F 2 /(l + R): 

d 2 a 2ita 2 ( , , y 2 \ „ , 

2 ( l ~y) + rr-s ^Q 2 ). (12) 



dxdQ 2 Q A x V 1 + R, 

In order to calculate the RC one has to use a realistic parametrization of the 
structure functions Fi(x,Q 2 ) over the full range of x and Q 2 . We have used several 
different parameterizations of the deep inelastic proton structure functions , including 
MRS D- and MRS DO parameterizations . 17 For region Q 2 < 5 GeV 2 the parameter- 
izations were taken at Q 2 = 5 GeV 2 and multiplied by factor 18 (1 — exp(— aQ 2 )) , 
a = 3.37 GeV~ 2 . In our calculations we have assumed R = 0. 

In the " Zeuthen analysis" the following x, Q 2 -bmning was used. There were four 
<5 2 -bins: 

Q 2 = 5 - 10, 10 - 20, 20 - 40, 40 - 80 GeV 2 (13) 
and four bins per rc-decade: 

log x = -4.0, -3.75, -3.50, -3.25, -1.0. (14) 



The kinematical region of the experimental measurement was defined by: 

157.5° <6 e < 172.5°, 0.05 < y e < 0.60 (15) 
for electron variables and 

157.5 <6 e < 172.5°, y h < 0.50, 9 jet > 10° (16) 
for mixed variables. 

Figs. 1,2 show the results for the radiative corrections in the terms of the mixed 
variables x m) Q 2 . which are calculated by the integrations over bins. The RC in the 
terms of the mixed variables, 5(x m ,Ql), are not small and can be reach 7 % at 
x = HT 2 . 

The "Zeuthen analysis" of the deep inelastic scattering has included also the 
analysis of the experimental data in the terms of the electron variables x e , Q\ corre- 
sponding to the binning ( 13), ( 14). Fig.3,4 show the results of the calculations of the 
RC in the bins x e , Q 2 e) that have been obtained with the program TERAD91. 

The corrections determined by the electron variables, 8(x e , Q 2 ), are known to be 
large at low x and dominated by hard photon emission from the incoming electron. 



From the comparison of Fig. 3 and Fig. 4 it is seen that the sensitivity of the RC to 
the choice of the structure functions is not small. 

The comparison of the results for the radiative corrections in the terms of the 
electron variables E' e) 9 e , that have been obtained 19 with event generator HERACLES 
( interface DJANGO 20 ) and analytical program TERAD93, is shown on Fig.5 for the 
bin 162.5° < 9 e < 167.5°. Here the variable u is defined by the equation 



and for the structure function F 2 the parametrization MRS D- was used and R=0. 
The results of the both programs are in good agreement. 

Now we will discuss the influence of the calorimetric measurement of the energy 
of the scattered electron on the size of the radiative corrections. 

Because the collinear radiation of photons in the final state of the reaction ( 5) 
is not resolved in the electromagnetic calorimeter, one must calculate the RC in terms 
of the calorimetric variables including the total energy of the electromagnetic cluster 



Fig. 6 shows the results for the radiative corrections including the effect of the 
calorimetric measurement of the final state photon radiation obtained 19 with the 
Monte Carlo HERACLES. For comparison the RC in the terms of the pure electron 
variables are shown also. 

In the program TERAD93 the possibility to calculate the RC in the terms of 
calorimetric variables, E vis ,6 e , has not been included yet, but with the help of the 
leading log calculations 21 one can estimate the size of the influence of the calorimetric 
measurement on the radiative corrections. 

The calorimetric measurement cannot distinguish between a single electron and 
an electron with accompanying photons. Therefore the final state radiation has exper- 
imentally no effect on the measurement of the kinematics of deep inelastic scattering 
and can be absorbed in the definition of the final electron energy of the process ( 1). 

This leads to the following way 19 to take into account the calorimetric mea- 
surement. We can calculate the contribution of the final state radiation in the leading 
logarithmic approximation and subtract from the total RC, that we obtain with the 
program TERAD93. 

The more detailed investigation of the calorimetry problem will be present in 
the separate publication. 
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Fig. 1. Radiative correction factor 8{x m ,Ql) (%)in (x,(5 2 )-bins for MRS D- parametriza- 
tion. 
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Fig. 2. Radiative correction factor 5(x m ,Ql) (%)in (x, Q 2 )-bins for MRS DO parametriza- 
tion. 
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Fig. 3. Radiative correction factor S(x e , Ql) (%)in (x, Q 2 )-bins for MRS D- parametrization. 





Fig. 4. Radiative correction factor 6(x e , Ql) (%)in (x, Q 2 )-bins for MRS DO parametrization. 
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Fig. 5. Radiative correction 5(u,9 e ) in ii-bins for 162.5° < 9 e < 167.5°. 
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Fig. 6. Effect of calorimetric measurement on size of radiative correction 5(u,9 e ) 
for 167.5° < 6 e < 172.5°. 



